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Abstract 

A detailed theory-versus-experiment comparison is worked out for Hj" intense laser dissocia- 
tion, based on angularly resolved photodissociation spectra recently recorded in H.Figger's group. 



' As opposite to other experimental setups, it is an electric discharge (and not an optical excita- 

r**- ■ tion) that prepares the molecular ion, with the advantage for the theoretical approach, to neglect 

without lost of accuracy, the otherwise important ionization-dissociation competition. Abel trans- 
formation relates the dissociation probability starting from a single ro-vibrational state, to the 
probability of observing a hydrogen atom at a given pixel of the detector plate. Some statistics 
on initial ro-vibrational distributions, together with a spatial averaging over laser focus area, lead 

a, 

to photofragments kinetic spectra, with well separated peaks attributed to single vibrational lev- 
els. An excellent theory-versus-experiment agreement is reached not only for the kinetic spectra, 
but also for the angular distributions of fragments originating from two different vibrational levels 
resulting into more or less alignment. Some characteristic features can be interpreted in terms of 
basic mechanisms such as bond softening or vibrational trapping. 
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I. INTRODUCTION 



The above threshold multiphoton ionization and dissociation of subjected to strong 
laser interaction, have revealed interesting nonlinear effects in angularly resolved kinetic en- 
ergy distributions of the photofragments, measured in experimental works covering the last 
decade Among these are the observations of very large increase (or some- 

times decrease) of the photodissociation rates originating from some vibrational states of 
the parent molecule at some specific laser intensities or, even more unexpectedly, misalign- 
ment effects in fragments angular distributions f|. The interpretation of such behaviors 
has been attempted by referring to some basic dynamical mechanisms evidenced through 
the light-induced adiabatic potentials describing the dressed states of the molecule-plus-field 
system. According to the frequency regimes, bond softening (in UV) or barrier sup- 

pression (in IR) [8( mechanisms tend to enhance the dissociation cross-section especially in 
the polarization direction of the laser. As opposite to them vibrational trapping (in UV) 
or dynamical dissociation quenching (in IR) [ijj], act as stabilization mechanisms, fa- 
voring misalignment in the fragments distributions. This complementarity has also been 
referred to, for laser control purposes of the chemical reactivity; namely by softening some 
bonds while hardening others Although very accurate quantum calculations in the 

frame of time dependent approaches have been carried out, with successful interpretations 
of dynamical behaviors in short-intense laser pulses, to the best of our knowledge, there is 
no a thorough and quantitative theory-versus-experiment comparison, up to date, the work 
of Kondorskiy et al. being a precursor in this direction. Basically two reasons can be 
invoked for the difficulty of such an attempt: only very few theoretical models take into ac- 
count the competition between ionization and dissociation processes leading, in very strong 
fields, to Coulomb explosions and only very few experimental works are conducted with a 
careful investigation of vibrational populations and sufficiently high momentum and angular 
resolution yielding accurate information about the dissociation of single vibrational levels. 

Experimental works on this system can be classified according to the preparation of the 
parent ion from the neutral molecule Ho. A first category collects experiments referring 
to optical ionization with a laser prepulse jjj, |2j, |3, |^ . The independence of the ionization 
and dissociation processes can not be experimentally controlled, and their competition is 
still an open question More recently, another kind of approach has been investigated 
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through ion beam experiments, where H^" ions are produced in a dc electric or plasma 
discharge that disentangle ionization and dissociation processes jl^i An accelerated 
and strongly collimated monochromatic beam is crossed at right angle by a focused 
intense laser beam. An advantage of the strong ion beam collimation is the reduction 
of the intensity volume effect; all ions being approximately irradiated by the same laser 
intensity (the validity of such approximation will however be discussed hereafter). Moreover, 
experiments conducted with low intensity pulses coupled to computational simulations of the 
resulting dissociation spectra, allow the determination of the population of the rovibrational 
levels of molecules in the beam. The neutral dissociation fragments (H atoms originating 
from photodissociation of Hj) are projected on a multichannel detector (MCD), whereas the 
charged particles (undissociated molecules and H + fragments) are extracted by deflection 
into a Faraday cup using an electric field. Excellent energy resolution (about 1%) allows 
the separation, in the circularly shaped patterns observed on the screen, the momentum 
projection of fragments almost originating from a single vibrational level . 

A model aiming in a quantitative theory- versus-experiment comparison, within the frame 
of the ion beam setup, has to fulfill the following requirements: 

i) the photodissociation process has to be accurately described in the center of mass 
frame by a wavepacket propagation under the effect of an intense radiative field, starting 
from a given rovibrational state. There is no need, however, to refer to any competition 
with ionization, as the experiment precisely disentangles these two fragmentation processes. 

ii) a geometrical transformation towards the MCD-plate has to be carried out, taking into 
account the macroscopic kinetics of the ion beam. This relates the total number of particles 
collected by a given pixel of the plate, during the whole experiment, to the previously 
calculated wavepacket, describing the evolution of an initial rovibrational state under the 
effect of a laser pulse of a given intensity. 

iii) although particular attention has been paid to the ion beam collimation in order to 
reduce the field intensity volume effects, a spatial average over the laser focusing area has 
to be carried, taking into account the different radiative couplings felt by molecules 
according to their geometrical position in the beam. This can be done through the use of 
some experimental measurements of the intensity distribution in the focus carried through 
a pinhole of 1 fim diameter 13j . 

iv) quantitative agreement also requires an averaging on the detector plate using some 
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windowing functions that simulate the resolution power of the detector. 

The organization of the paper follows these achievements in Section |H] The results and 
their interpretation are presented in Section ITTT1 with a thorough discussion of the role of the 
intensity volume effect. An excellent theory- versus-experiment agreement is obtained not 
only for the kinetic but also on the angular distributions of the photofragments. Section IV 
is devoted to some conclusions and perspectives. 

II. THEORY 

Referring only to two radiatively coupled Born-Oppenheimer electronic states; namely the 
ground (lsa 9 ) and the first excited (2pa u ), an accurate wave-packet propagation method using 
the split operator technique is described in detail in ref . [3| . For the sake of completeness, 
we give hereafter a brief summary of the method, introducing the corresponding coordinates, 
operators and quantum numbers. The emphasis is rather put on the way to relate the 
quantum information content of the wavepacket to the observed momentum projections of 
the neutral photofragments H resulting from a rovibrational distribution of parent ions 
excited by a laser source of given spatial distribution. 

A. The wavepacket propagation 

In the laboratory frame and using spherical coordinates, the total molecule-plus-field 
Hamiltonian is written in terms of a two-by-two operator matrix: 



R is the diatomic internuclear vector. R, 6 and <p designate the internuclear distance, polar 
and azimuthal angles of R with respect to the laser polarization vector e, respectively. As 
is usually done, a functional change on the wavepacket: 



H{R, 6, 0; t) = T R + T e + T + V(t). 



(1) 



1 



(2) 



R 



aiming in a simplification of the radial part of the kinetic operators, leads to: 



T R = -1 



1 d 2 
2M dR 2 



(3a) 



4 



2MB 2 sin 2 9 d(f) 2 



with 1 the identity (2x2) operator matrix. Atomic units (^=1) are used in Eqs(JH} where 
Ai designates the reduced mass. The time dependence arises in the non-diagonal terms of 
the potential energy operator matrix V through the radiative couplings: 



where fi(R) is the transition dipole moment and £ (t) is the laser electric field amplitude, 
given as product of a pulse shape e(t) times an oscillatory term involving the carrier wave 
frequency to: 



Note that the cos 9 in Eq. (J1J results from the dot product of the transition dipole vector 
(parallel to R) times the laser polarization vector e. 

The diagonal elements V\(R) and V^-R) of V are nothing but the BO curves of the ground 
(label 1) and first excited (label 2) states of H^~. V\, V 2 and /i are obtained in the frame of the 
Born-Oppenheimer approximation, at the zero order level with respect to the ratio m e /m 
of the electron to the proton masses. Using spheroidal coordinates, it is well known that 
the Schrodinger equation can be written as two eigenvalue equations Q , which have 
been numerically solved here using the shooting method j^jj. The potential energy curves 
have been computed in the range < R < 200 a.u., with a numerical accuracy checked to 
be better than 10 -12 a.u. The mass ratio m/m e has been taken as m/m e = 1836.152701. 
Finally the dipole matrix element \x between the lsa g and 2pa u states has been obtained by 
numerical integration of the wave functions, at the same level of numerical accuracy. 

The time dependent Schrodinger equation (TDSE) describing the wavepacket propagation 

is: 



V 12 (R,9,t) = n(R)£(t) cos#, 



(4) 



£{t) = e(t) cos out. 



(5) 



i^*(R, 9, 4>; t) = H(R, 9, 0; t)*(R, 9, 0; t) 



(6) 



with, as an initial condition: 





reflecting the fact that at time t = 0, only the ro- vibrational levels of the ground electronic 
state are populated. The eigenfunction $! precisely corresponds to such a state with quan- 
tum numbers g, v, N, Mm (electronic ground, vibrational, total and e-projected rotational) 
and is given by: 

$i(i2,0,0;O) = Xg,vAR)PN lf (<*xO)e iMN *. (8) 

P n In (cos8) is the (N,Mm) Legendre polynomial, whereas the radial part is defined as the 
solution of the time-independent Schrodinger equation: 

1 



d2 + 



E 



Xg,vA E ) = °' 



(9) 



2MdR 2 v y 2MB 2 
The motion associated with the azimuthal angle <fi remains separated under the action of 
the 0-independent V, such that is a good quantum number describing the invariance 
through rotation about e. 

The propagation u sing the split-operator technique has been described in full detail in 



previous works 



ralia 



221 ] . The peculiarity of odd-charged homonuclear ions is their linearly 



increasing dipole moment with R, leading to asymptotically divergent radiative couplings. 
We take them into account by splitting the wavefunction into two regions, an internal and 
an asymptotic one. The latter is analyzed by a generalization of the Volkov type solutions 
23], while the numerical propagation on the former is performed by Fourier transform 

n n 

methodology [24| with the implementation of a unitary Cayley scheme for Tq |22j . 



B. From wavepacket to observed spectra 

The main concern of this paragraph is to relate the experimental observable, i.e. the 
probability distribution of hydrogen atoms resulting from H^" photodissociation, as recorded 
on the multichannel detector (MCD), to the asymptotic part of the wavepacket $(R, 6, <p] t) 
solution of Eq.©. By asymptotic we mean large internuclear distances R for which the 
molecule is considered as dissociated without the possibility of a recombination process. 
To the best of our knowledge such a correlation has not rigorously been attempted in the 
literature. So far, the interpretation of general tendencies of photodissociation spectra re- 
ferring to basic mechanisms, has rather been conducted by angularly resolved kinetic energy 
distribution given by: 

V(k, 6, <p) = lim $(k, 6, <p;t) 2 , (10) 

t— >oo 



where 

1 r°° 

$(k, 9, 0; t) = —= \ ®(R, 9, 0; t)e~ ikR dR. (11) 

V27T J-oo 

is the Fourier transform of $ over the scalar variable R, (i.e. not over R taken as a vector). 
The argument retained by doing so, is that asymptotically, due to R~ l type of behavior 
in the kinetic operators Eqs (j3b|3cj) . angular dynamics is not affected at large internuclear 
distances. Note that in this paragraph, for the sake of simplicity, we drop the labels of $ 
depicting initial state quantum numbers (v, N, M^). 




FIG. 1: The photo-dissociation experiment through the H2 photoionization 

To reach a comparative level of understanding, we are now describing the two families of 
experiments. Pertaining to the first family are photodissociation experiments where both 
photodissociation and photoionization steps are laser induced Starting from 

neutral H2, in its ground electronic and vibrationless state X(t> =0), a multiphoton excitation 
leads, through the EF intermediate excited electronic state, to the H^" ground state with a 
distribution of rovibrational levels. Dissociation follows the absorption of additional photons 
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and is very fast as compared to the relative motion of the parent ion in the laboratory 
frame. Whence the photofragments are well separated, H + ions are extracted (accelerated) 
through an electric field and collected on the MCD plate. A schematic view is provided in 
figure d Photodissociation occurs, as a fast process, at the origin O of the laboratory frame, 
at a time which is taken as t = 0. The laser polarization vector is along the z-direction, ri 
and T2 are the vectors pointing H and H + . A further step is the extraction of the proton 
H + by an electric field applied along the y-direction towards the MCD plate positioned at 
a distance 00'=D from the origin. The detection occurs on a pixel M defined by its polar 
coordinates (p, a) on the MCD surface (or by r with respect to O) that H + is reaching 
after a time of flight t, with velocity v. It is worthwhile noting that this last step is just a 
mapping of the photofragment onto the detector (without dissociation during time t). The 
vector transformation relating the proton H + position (R, 6, </>) in the center of mass frame 
to the pixel M (p, a) on the detector is known as the Abel transformation j^ . 




FIG. 2: The photodissociation experiment based on the ionization of H2 using discharge source 



S 



A different situation prevails in the ex peri ments of the second family where an electric or a 
plasma discharge ionizes H 2 into [l3lll(| . The resulting ion beam is strongly accelerated 
by an electric field and is crossed at t = by the laser beam at a point O of the laboratory 
frame. The description of such experiments, as illustrated in figure has to combine two 
motions; namely, the translation of the center of mass G in the laboratory frame along u v 
(unit vector along y) with velocity v and the nuclear separation (dissociation) in the center 
of mass frame. The hydrogen atom H resulting from photofragmentation is collected at the 
pixel M of the detector. It is to be noted that M is positioned with respect to the laboratory 
frame with a vector r, corresponding to r\ at time t, when H reaches M. 

As, our concern is the quantitative interpretation of photodissociation _spectra obtained in 
H.Figger's group using an electric discharge to induce ionization [3,0,01; emphasis is put 
in the following on a thorough description of the kinematics of the second family experiments. 
The quantity which is measured, is nothing but the number of hydrogen atoms dN collected 
on each pixel M (p, a) at infinite time. This can ultimately been related with the time 
integral of the flux of the current density j(p, a, t) of H orthogonal to the area dS = pdpda 
of the finite size pixel M, as 

POO 

dN(p,a) = NdS j(p,a,t) -Uydt, (12) 
Jo 

where N is the total number of H photofragments. The flux in Eq.(|12j). involves an averaging 
over the positions of all protons H + that are not detected in the experiment jl^ |: 

j(t) = - [dr 2 lm[**(R,RQ-,t)V n *(R,Ro;t)). (13) 

where ^(R,Ra;t) is the overwhole wavepacket describing the combined molecular inter- 
nal i2-motion and center of mass i2o-motion. Im stands for the imaginary part. Frame 
transformations defining Rq and some vector relations directly related with figure El are 
gathered in an Appendix. Separation of the photofragments relative motion described by 
^(R; t), from the motion of the center of mass described by &g{Rg] t), leads to the following 
representation of the total wavefunction: 

^(R.R G :t) = ^(R:tH> ( AR G :n. (14) 

Using the frame transformations Eqs (jA-ipA-2)l together with Eq.(J2J) one has: 

= *^Z^ 4o (f!l±5); t ). (15) 



Concerning the calculation of the gradient V ri involved in Eq.([13jh we note that the flux 
has to be evaluated at large R with as consequences: 

V r MR; t) = V R $(R; t) ~ u R ^{R- 1) (16) 

with Ur the unit vector along R (Eq. (|A-6|) ) and 

V ri $ G (i2c; t) = ^V^aiRa; t) (17) 

The approximation involved in Eq. ljlfjjl results from the neglect of all angular derivations 
due to their occurrence with coefficients decreasing faster than We proceed now to a 
quasiclassical approximation for the description of the center of mass translational motion, 
with two implications: 

(i) Rq — vUyt has a corresponding wavevector Kg — mvu y and the application of mo- 
mentum operator — iVR G to $g(-Rg) simply result into 77w$g(.Rg). When this is done at 
the level of Eq(JT3J) one gets: 

d 



If 1 

j(t) = — / dr 2] rrlm 

m J \T\ — Ti\ 



§*(R] t)^-$(R; t) + imv\$(R; t)\ 2 



\$ G (RG;t)\ 2 . (18) 



(ii) No wavepacket spreading is allowed for $g(.Rg;£) which is localized with an envelope 
behaving as a 5-like function, i.e.: 

\<5> G (RG)\ 2 -S(RG-vt). (19) 

The integration over r 2 (with dr2 = 2dRo) finally leads to: 

(20) 



mR z 



$*(R; t)^(R; t)u R + imv\${R] t)\ 2 Uy 
oK 



R=2r\ —2vtu y 

with a rather intuitive interpretation of the two components of the flux. The first i.e. 
$*(i2; t)-Jj^&(R; t)u R is merely the current density generated by the expanding wavepacket 
in the center of mass frame, whereas the second corresponds to the current associated with 
a density travelling with a velocity v along u y . The calculation can be further conducted 
analytically by deriving an asymptotic (i.e. R —> oo, t — > oo) expression for $ [15]. This is 
done using a time evolution expression involving the Fourier transform Eq. ljll)) . Actually, 
one has for large R, where the potentials can be considered as constant and after the laser 
is turned off: 

6, <p; t) = e- ikH/m ^(k, 6, <P) (21) 
10 



the solution being induced only by the radial part of the kinetic energy. Returning back to 
the wavepacket in the coordinate space: 

1 



dk$(k, 9, (f))e- lk * /r V fcH 



'2-K j-oc 

and replacing $ by its expression Eq. (fTT| . one gets: 

dR^(R',9,(l))e im{R - R ' )2/4t . 



m \V2 r°° 
Jo 



Ai-nt 



Expanding the independent part of the exponential as: 



im / p d'\2 

g 4t ( n n > = g 4t g 



imR? imRR 1 



and observing 



g At 



that for large t: 



lim 

t— >oo 



imR 

e 4t 



1 



0. 



an asymptotic expression is obtained for $(-R, 0; £): 

/ m \ V 2 imR 2 - / mR 



V2i£/ 



2f 



0. 



(22) 



(23) 



(24) 



(25) 



(26) 



While recasting Eq.(|26J) into Eq. (j2U|) . a rather simple expression results for the asymptotic 
flux: 



m 



m 



2Rt 2 



2vt 

UR + — Uy 



2t ' ' 



(27) 



H=2(ri— vtUy) 

The calculation of the projection of j on M y (cf Eq ll2j) requires the vector relation of 
Eq. ()A-7|) that finally leads to: 



mD 



1 



j u b 



mR 
~2T 



(28) 



* 4t 2 [p2 + (£> - ut)2]2 

where i? depends on t as given by Eq. (|A-3|) . 

The final step is to transform the time integration of j ■ u y , involved in Eq. (jl2j) . into an 
integration over the kinetic momentum. We proceed to a change of variable: 



r 9 / ^ x9 n 1 /9™ R TTIV 

k=[p 2 + (D-vt) 2 }^ 2 — = - — , 



(29) 



D 2D 

the physical meaning of which will be clarified hereafter. Straightforward calculations show 
that Eq.(|29jl can be inverted as: 

forfe[0,f] (fee^^ + mV)^]) 
(fc2 ' fcg)1/ ^ for * 6 [f , +oo] (A; G [*„, +00]) 



#(1 



D 



1 + 



(30) 
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upon the introduction of the notation 



k p = mv— 



and leads to: 



dt = =f 



D 



k 



:dk. 



(31) 



(32) 



mv 2 ^2 _ /^j 1 / 2 

The =F signs correspond to the two time intervals depicted in Eqs (j30j) . The time dependent 
argument of $ in Eq. ()26|) can than be expressed using the two variables k fEq. ()29|) ) and k p 
(Eq.(jSH)) as: 



mR 
~2A~ 



m 

T 



V + (D- vt) 



° k=k f 1T M)l /2 ' 

■yt \ mv 



(33) 



The experimental conditions, are such that the velocity v of the molecular beam is much 
greater than the fragments relative velocity. We can thus consider {k 2 — k 2 ) 1 / 2 /mv as 
negligible when compared to 1, taking into account that D is much larger than p. The 
resulting approximation, namely: 

mR 



t 



~ — and 

v 



2t 



k 



(34) 



merely means that the time needed for a fragment to reach the pixel M(p, a) is approximately 
the same as the one needed for the center of mass G to reach the center O' of the detector. 
In the framework of this approximation, the meaning of k p = mpv/D ~ mp/t (defined by 
Eq. (}3*Tj) ) is also clear: i.e. the projection of the kinetic momentum k on the detector plane. 
Finaly we obtain for the time integrated flux: 

2 



j ■ Uydt 



1 1 

m v 
AD 2 



+ / 

kp J kp 



$(fc,M) 



k (k 2 - k 2 ) 



1/2 



dk 



(35a) 



2 2 

m v 
2D 2 



M 2 



dk. 



(35b) 



'kp k(k 2 -k 2 ) 

where the upper bond of the second integral in Eq. (j35a|) has been extended up to +oo 
considering that \$>(k, 9, (j))\ = for k > mv, which is equivalent to state that the center of 
mass kinetic momentum 2mv is much larger than the relative momentum of photofragments 
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k. Recasting Eqs (j35j) in Eq. (fT2|) . taking into account cylindrical symmetry over and 
calculating the pre-integral factor as: 



m 2 v 2 mvp mvdp 

-^-dS = — — da = kpdkpda. 



one finally gets: 



dN(k p , a) = N kpdkpda - 



k (k 2 - k 2 ) 



1/2 



dk. 



(36) 



(37) 



The dependence over 9 of the right-hand-side of Eq. (}3Tj) has to be expressed in terms of a, 
referring to the frame transformation Eq. (JA-9J) 

k, 



cos 6 = cos a 
k 



(38) 



in such a way that, ultimately dA" is written only in terms of the variables k p and a, with 
the parameters v and D characterizing the experimental setup: 

2 

-dA;. (39) 



dN(k p , a) = Nk p dkpda- 



$ (k, axccos(k p / k cos a)) 



1/2 



k (k 2 ,. j 

The probability to record a hydrogen atom on the surface element dS (pixel M) located at 
p, a on the MCD (with a kinetic momentum k p ) is obtained by a proper normalization: 

2 

-dk. (40) 



P(k p ,a)dS=^dN(k p ,a) = ^- 



$ (k, aiccos(k p / k cos a)) 



k [k 2 - kfj 



1/2 



It is interesting to note that the two probabilities V(k, 6) (given in Eq lTT7j) and P(k p , a) 
(Eq. p0|l ) are simply connected by: 



P(k p ,a) 



V(k, 



k (k 2 - k 2 p ) 



1/2 



dk. 



(41) 



Two remarks are in order: 

i) both equations Eq. (j4*Uj) and Eq.(jHJ involve a singularity at k — k p . This difficulty can 
be overcame by a partial integration leading to: 



P(k p , a) = —— arccos 
2 lip 



V(k,9) 



k=oo 



arccos 



fUw)d* (42) 



k=k P 2k p Jkp 

The integrated term in the right-hand-side of Eq. (j42j) is null, due to the fact that 



k=oo 



0. As for the total derivative with respect to k of V(k,9), it results into: 
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jL-p(k 0) = —+ fcpcosa &P 

dk 1 ' } dk k(k 2 - k 2 p cos 2 ay/ 2 06 ' 



When recasting Eq.()43|) into Eq. 

P(k p ,a) = 



one obtains: 



arccos 



p Jk t 



dV 



k p cos a 



k ) \_dk k(k 2 — k 2 cos 2 a) 1 / 2 89 



dk 



(43) 



(44) 



For k ~ k p , the singularity in the coefficient of S may only occur for a = 0. This is 



fortunately compensated by the aiccos(k p /k) term of Eq. (j42j) . Actually expanding Eq. 
in terms of powers of (1 — k p /k) one ends up with a non-singular behavior, i.e. ^ or ^ e 

integrand in the vicinity of k p . 

ii) Despite the fact that the experimental situation we are describing is not amenable to 
a simple mapping on the detector plate of a photodissociation that had already occured in 
the center of mass frame (as a figure Q), it turns out that Eq. (j4T?|) can finally be recast in 
terms of the commonly used Abel transformation |2a |: 

2" 



P(k p ,ot) = -A 



$ (k, arccos(/c p //c cos a)) 



where A is defined as Ml 



A[f{k)]{x) 



k 



f{k)k 



(45) 



:dk. 



(46) 



fx {k 2 -x 2 ) 1/2 

In connection with this, it is worthwhile considering the full Fourier transform of the wave 
function describing the relative motion (in contrast to the one carried in Eq lll|) : 



(2tt) 3 / 2 




dffl(R)e 



ikR 



(47) 



HP 



It can be shown by following the same derivations as Eqs (j21H27)|) jl5ll27|. that asymptotically 
(i.e. for t — > +oo and R — > +oo), one has: 



which, combined with Eqsfl2*landl2*6^). implies that: 

$(M 12 



(48) 



(49) 



The probability in Eq. f4*5|) appears now as the Abel transform of \l/ (k, arccos(fc p //ccosa)) 
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C. Rovibrationally averaged spectra 



The probability distributions which are calculated in the previous paragraph refer to 
a given initial state (g, v, N, Mjy) involved in the determination of $(t = 0) through 
Eq.(jBJ) such that the quantity resulting from Eq.fJSJ) is actually P v ,N,MN(k P ,o:) using a full 
notation. An averaging over the initial ro-vibrational populations is thus required to reach 
the experimental spectra. As the rotational states N are (2N + 1) times degenerated, a 
summation can be carried out over M N , leading to: 

1 N 

P v ,N(k p ,a) = 2N + 1 c M N Pv,N,M N {k P ,a), (50) 

M N =0 

where cq — 1 and cm n = 2 (for Mn ^ 0), due to the fact that the total Hamiltonian does 
not depend upon the sign of Mjy. The homonuclear diatomic character of H^" implies a 
total wavefunction (accounting for the nuclear spin) that is antisymmetric with respect to 
the interchange of identical nuclei. To ensure such a property the total nuclear spin number 
T must be either (associated with even N), or 1 (associated with odd N). Due to very 
rare singlet (T = 0) - triplet (T = 1) transitions, molecular hydrogen mainly consists of two 
distinct species: parahydrogen (T = 0) and orthohydrogen (T = 1). The occurrence of ortho 
states is three times more probable than the para ones |l7[. This nuclear spin statistics is 
accounted for, by a weighting coefficient g^, affecting the rotational populations N, such 
that 

{1 for even iV 
(51) 
3 for odd N 

Apart from this factor, rotational populations are also thermally weighted, according to a 
Boltzman distribution described by a rotational temperature T v depending on the initial 
vibrational level. The weighting coefficient is given by: 

T AE(v,N)] 

b N = ex V — . (52) 

where hp stands for the Boltzman constant and AE(v, N) for the rotational energies resulting 
from the solution of Eq. : 

AE(v, N) — E(v, N) — E{v, 0). (53) 
The rotationally averaged probabilities resulting from these considerations are 

P v ( k p , a) = — ^2 bNgNPv,N(k p , a) , (54) 

N 
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where Q v is a normalization factor: 

Qv = ^2 h N9N- (55) 

N 

The comparison with experimental spectra has also to take into account initial vibrational 
populations a(v) (i.e. as they result from the electric discharge acting over H 2 , prior to the 
laser excitation): 

P(k p , a) = — ^ a(v)P v (k p , a), (56) 

V 

with 

Q = Y j a(v). (57) 

V 

We note that the informations concerning the initial vibrational distribution a(v) as well as 
the corresponding rotational temperature T v , have to be provided by experimental measure- 
ments. 



D. Laser spatial intensity averaging 

Although particular attention is devoted in the experimental measurements for obtaining 
a well focussed ion beam, the molecules are actually excited by different field amplitudes 
according to their positions (x, y), due to a non-homogeneous spatial intensity distribution 
I(x, y) in the laser beam. It turns out that an average over these non-homogeneities has 
a basic importance when attempting a quantitative interpretation of experimental data, as 
will be clear in the next section. The average implies a double spatial integration over the 
variables x and y (see figure [BJ): 

/l r+co 
dx / P(k pj a;I(x,y))dy. (58) 
L J — oo 

x being limited to a finite interval 2L measuring the width of the ion beam. The integrand 
itself is nothing but the probability calculated in Eq. ()56j) with as an additional argument the 
intensity I(x,y), for which this probability is calculated, i.e. P(k p ,a;I(x,y)). For parity 
reasons Eq. (j58|) may be also written as: 

/•L r+oo 

P(k p ,a) = A dx P(k p ,a;I(x,y))dy. (59) 
Jo Jo 
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A gaussian shape is assumed for the 2D behavior of the laser within its focus area: 



I(x, y) = I exp 



x 2 ' 






^ x _ 


exp 


- r y- 



(60) 



where r x and r y are the radii of the focal area, such that I(r x , r y ) = Io/e 2 . These parameters 
are obtained from the experimental setup |l3j] as: 

A/ 



2nb 



(61) 



x,y 



where A is the laser wavelength and / is the focal length of the parabolic mirror focusing 
the laser beam. b x and b y correspond to the extensions in each direction x and y where 50% 
of the energy is dissipated. The peak intensity is calculated whence the total pulse energy 
Eq and an autocorrelation time t ac have been measured: 



2V2 ln2£ 

Tl^/ 2 T x Tyt ac 



(62) 



The double integration involved in Eq. (j59jl can be conducted, by integrating first over y, 
referring to a variable change: 



y = r y 



dy = - JL 
y 2 



-In 



-hi — 



1/2 



where I x = io ex P [~( x / r x) 2 }- The result is: 

P(k p , a; I x ) = 2r y 



u P{k p , a; I) 



1/2 



dl. 



(63a) 



(63b) 



(64) 



Two singularities affect such an expression; namely at I = and at I = I x . The first has no 
consequence, as for 1 = 0, P(k p , a; I = 0) = 0. The second can be avoided when integrating 
by part: 



P(k p , a; I x ) = 2r y 



-2^TJT x P(k p , a; I) + J X 2^T/T x ^jP(k p , a; I)dl 



(65) 



An identical procedure is then applied for the second integral over x. The final result is 



P{k p , a) = 2r x 



Ix=Io 



-2y/I x /I P(k p ,a;I x 
where II is the field intensity at (x = L, y = 0) 
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d 



+ / 2y/I x /I —P(k p} a;I x )dI x 
Ix=Il Ji l dlx 



(66) 



III. RESULTS 



This section presents the results of the simulation and interpretation of experimental data. 
Among the experimental results obtained in H.Figger's group ^ll^El three are retained. 
The laser intensities Iq have been slightly adjusted with respect to laboratory measurements 
of the total pulse energy E and autocorrelation time t ac , to fit experimental spectra. The 
resulting parameters are collected in Table |U The laser pulse carrier wavelength is A=785 
nm. In the calculations the intensity spatio-temporal distribution is assumed to be: 



E (mJ) 


tac (fs) 


I (TW/cm 2 ) 


0.3 


228 


7.5 


0.5 


234 


9.5 


0.7 


240 


16.0 



TABLE I: Total pulse energy Eo and autocorrelation time t ac of the laser field. Iq is the maximal field 
intensity value. 



I(x,y;t) = J exp 



x 2 








" 2t 2 ' 




exp 




exp 




^ x _ 




- r y- 




. w t_ 



(67) 



with the relations between the parameters as given by Eqs()611 IB^j) . The widht of the 
molecular jet is L=50 /im, its velocity is t>=10 6 m/s and focal area parameters values are 
b x =2.6 mm, b y =2A mm, /=1000 mm, resulting into r x =48fim, r y =52/im. In the calculations 
described below the parameter w t = t ac /2\ / ln2, which define the laser pulse temporal shape, 
have been taken equal to 140 fs. 

Figure El displays three-dimensional representations of the dissociation probabilities as a 
function of their angular (a) and kinetic (k p ) distributions. The upper diagram corresponds 
to the experimental results P e [25] for the laser excitation parameters indicated on the last 
row of Table |U The lowest two diagrams give the calculated spectrum (P c ) and the absolute 
value of the difference \P e — P c \, for the same laser characteristics, at the same scale. The 



normalization is such that: 



7T/2 



P e ' c (k p , a)da 



(68) 
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The successive peaks that are obtained correspond to photofragments arising from different 
vibrational levels v of the parent ion H^~. The energies of the levels are positioned in figure 
0]on the rotationless dressed molecular potentials resulting from the diagonalization of the 
radiative interaction at fixed angle 8 = 0. 

The most important peak (at a = 0) corresponds to v=7 and is followed in decreasing 
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1.5- 




FIG. 3: (a) - The three dimensional representation of the experimental result corresponding to Eo=0.7 
mJ (see the third line in Tab. {I}; (b) - The corresponding calculation result; (c) - The difference between 
experimental and calculated spectra. 

order by the peaks assigned to w=8,9,10. The peak resulting from the dissociation of v =6, 
with a much smaller contribution is hidden by the peak v—7, whereas those resulting from 
v =11,12 are in the blue-wing of v=10. It is interesting to note that, on the lowest diagram, 
the largest error affects the peak resulting from v=9, all others being well represented. This 
is to be relationed with the particular energy of v=9 (see figure HJ) very close to the avoided 
crossing of the dressed potentials. The characteristics of this region being very sensitive to 
the laser spatial and temporal intensity distributions, even small deviations with respect to 
experimental evaluations may lead to appreciable differences explaining figure Efc). 

In figure El we show the four main steps to obtain the photofragments distribution P c , 
which may be compared with the experimental one. On each step we plot the cut of the 
resulting distribution at 9 = for V(k,6), or a = for P(k p ,a). The upper panel gives 
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FIG. 4: 2D view of the dressed adiabatic potential curves of Hj" (solid lines) (with a continuous wave laser 
of wavelength A=785 nm and intensity 7o=l-6xl0 13 W/cm 2 ) and the corresponding diabatic BO electronic 
states (dashed lines). Are also indicated the eigenvalues of the field- free vibrational levels. 

the photodissociation probabilities starting from individual vibrational levels of Fhj", as cal- 
culated in the molecular frame for #=0 and as a function of k, for laser characteristics 
corresponding to the last row of Tabled The vertical lines illustrate the theoretical energies 
of the vibrational levels of in a field-free situation. As expected, from the examination of 
figure IH the maxima of the peaks with v < 9 are shifted down and these of v > 9 are shifted 
up, due to the radiative coupling. As to the height of the successive peaks, a decrease from 
v=Q to v =9 is observed. This, however, does not mean that v—9 is less dissociated than 
v =6, as the information which is displayed concerns a cut at angle 6 = 0. 
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The major effect, when attempting a theory- versus-experiment comparison, is the role 
played by the spatial intensity distribution of the laser, that so far, has been neglecte d by 
referring to the large laser focal area with respect to the diameter of the ion beam jl2|. 
Spatial averaging brings into interplay molecules interacting with radiative fields having 
intensities less than the maximum value Iq. This may lead to very large effects on some 
vibrational levels as is seen in figureElpanel (b). More precisely, the peak associated with v=6 
is nearly washed out, whereas those describing photodissociation starting from v=9 and v =11 
seem to be enhanced. Only the laser maximum intensity leads to a barrier lowering (bond 
softening) mechanism for v=6. When a spatial average is carried out, with the inclusion 
of lower intensities, the photodissociation from v =6 is severely inhibited due to very low 
tunneling, which explains the flat behavior of V % v {k) for t>=6.The vibrational states v =7,8 
are also affected by this effect but in a lesser extend as they are closer to the top of the lower 
adiabatic potential barrier. The level v=9 is again in a particular situation, energetically 
lying on the very top of the barrier. In other words, its photodissociation is not inhibited by 
any spatial intensity averaging, and this is why it leads to a narrower and higher peak, than 
the ones originating from v =7,8. The narrowing of the peak, in particular, is in relation with 
the fact that only a limited energy range corresponds to efficient dissociation within the open 
gate between the lower and upper adiabatic potentials that gets narrower with decreasing 
intensities. The behavior of v—11 deserves particular interest, as its photodissociation is 
rather through a vibrational trapping mechanism involving the upper adiabatic potential. 
Lower the field intensity and lesser is the efficiency of this trapping. The spatial intensity 
averaging of the laser leads as a consequence to better relative dissociation from v—11 
resulting into a high peak. 

Figure El panel (c) displays the dissociation probabilities P^(k p ) as functions of k p after 
the Abel transformation Eq.([4ip. The basic difference with V l v {k) (panel b) is the rise of 
long range red tails of the individual peaks, especially for v > 9. This is due to the nonlinear 
features of the Abel transformation, mixing up a whole range of ^-dependent probabilities 
for a single a. Less aligned fragment distributions resulting from v=9,10,ll present tails 
that are much more marked than the ones coming from v—7 and 8. 

The following step for building the experimental observable is the convolution with the 
detector resolution window which is taken as a square gate of 0.07 a.u. in kinetic momentum 
units, corresponding approximately to a pixel size of 70 [im. This as expected, results into 
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FIG. 5: Successive steps for a theory-versus-experiment comparison of photodissociation spectrum for a 
laser energy 0.7 mJ (last row in the Table©. Panel (a) gives the individual probabilities of initial vibrational 
levels v—6,...,12 in the molecular frame, for 9 = 0, and as a function of k. Panel (b) takes into account the 
spatial intensity distribution of the laser. Panel (c) displays the intensity averaged probabilities after Abel 
transformation bringing them in the laboratory frame, for a = and as a function of k p . Panel (d), same as 
c but after convolution by the detector resolution window. Panel (e) sums up all individual v contributions 
and compares with the experimental spectrum. 
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the smoothing of very sharp structures such as the peaks associated with v=9 and 11 (pannel 
d). 

Finally the theoretical spectrum is obtained as a sum of partial vibrational distributions 
with weights corresponding to the vibrational populations as given by Eqs (j56Hf)Tj) . Panel (e) 
of figure displays the resulting kinetic energy spectrum which is directly compared with 
the experimental one. An excellent agreement is obtained, the most noticeable differences 
occurring in the vicinity of v=9, which corresponds to an energy region particularly sensitive 
to possible inaccuracies related with the spatial intensity averaging of the laser. 

The information we get from figure can be summarized as follows: apart from the 
geometrical Abel transformation which is needed to bridge the dissociation probabilities 
evaluated in the ik-6) frame, to the photodissociation spectra as recorded on the detector 
plate, one has to take into consideration basically two additional facts, when attempting a 
quantitative interpretation of the experimental data: 




FIG. 6: The laser field intensity distribution over the molecular beam. 

i) The first is the spatial intensity distribution of the laser. Figure |B1 displays a three- 
dimensional view of the relative spatial extensions of the laser focal area and that of the 
molecular beam as is actually the case in the experiments. Clearly, all the molecules are not 
subjected to the same intensity at a given time, requiring thus a spatial averaging, the role of 
which is one of the most striking. Figure[7|gathers the spectra on the detector plate, for a = 
and as a function of k p for two models: namely, with and without the spatial averaging over 
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FIG. 7: The dissociation probabilities calculated with (dashed line) and without (dotted line) averaging on 
the laser intensity distribution, versus experimental data (solid line). 

the laser intensity distribution. The results are compared to the experimentally recorded 
data. A huge decrease affects the spectrum in the momentum region extending from k p ~ 
3 to 5 a.u. when averaging over the field intensities. This precisely corresponds to the 
contributions of vibrational levels v =6,7 well protected against potential barriers that are 
high for lower intensities taking part in the averaging process. Thus the spatial averaging 
turns out to be crucial when comparing with experimental spectra. 

ii) The second is an accurate knowledge of the field-free vibrational populations of the 
parent ion Hjj", which take part in Eq. (|5l)|) through the function a(v). Figure |H1 displays in 
terms of histograms the relative vibrational populations of levels v =6,..., 12 as they result 
from an estimation based on similar discharge experiments (^J. They are actually subjected 
to errors presumably within 10 to 15% in relative values. A calculation based on them leads 
to the spectrum illustrated in figure |H1 which basically disagrees with the experimental one 
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FIG. 8: Vibrational levels populations as estimated by Sandig in 



and fitted in the present work. 



over a region close to k p ~ 5 to 6 a.u., corresponding to the most important peak. However, a 
nice agreement is obtained after some small modifications of the vibrational populations, not 
exceeding reasonable error bar limits and remaining within the overall decreasing behavior 
for high f's, as is plotted in figure El It is worthwhile noting that this adjustment is done 
once for all, for given laser parameters (0.7 mJ of total energy) and is used hereafter for all 
other theory-versus-experiment comparisons. 

The spectra are gathered within the frame of two one-dimensional representations: either 
as a function of the kinetic momentum k p or as a function of the angle a on the detector. 
Figure CHI gives the cuts (at a=0) as a function of k p , for three laser fields, whose character- 
istics are precisely the ones indicated on Table HJ It might be noted here that peak intensity 
value Jo calculated using Eq. (jH2J) strongly depends from the experimentaly measured param- 
eters b x , by and E Q . In third column of Table H] we give the intensities corresponding to the 
best agreement between experimental and calculated spectra presented on figure [TUJ This 
adjustment is necessary to reproduce correctly the left part of the spectra, which is partic- 
ularly sensitive to the laser intensity, as can be seen in figure We emphasized that this 
adjustment has only been performed for one-dimensional spectra corresponding to (a=0). 
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FIG. 9: Dissociation probabilities calculated using Sandi. 
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brational levels populations (dotted line) 



and the modified ones (dashed line) as compared with the experimental data (solid line) . 

For the same pulse duration the laser intensities are ranging from low to medium and strong, 
following the panels a,b, and c. Three features can be emphasized: 

i) Three major peaks are obtained, corresponding to the dissociation involving v—7, 8 
and 9 levels, positioned in this order in the region k p ~ 5 to 8 a.u. The theory-experiment 
agreement is good not only for the positions but also for the relative heights of these peaks; 
the most noticeable difference affecting again v =9, more sensitive to an accurate evaluation 
of the spatial intensity distribution; 

ii) The strongest field (E =0.7 mJ, panel c) reveals the rise of an additional peak at the 
position of v =6. This is related with the bond softening mechanism, where the radiative 
coupling induces an important adiabatic barrier lowering, large enough for the population 
initially in the bound state v=6 to escape through tunneling. Excellent agreement with 
experimental results is obtained for this bump in the spectrum ; 

iii) The blue tail of the spectrum extending above k p ~ 8 a.u. corresponds to the photodis- 
sociation of initial populations on v =10,11 and 12, which due to possible vibrational trapping 
effects, is less efficient. Here again excellent theory-experiment agreement is reached. 
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FIG. 11: Angular distributions of photofragments for: (a) E o =0.3 mJ; (b) E =0.5 mJ; (c) E =0.7 mJ. 
Solid lines correspond to experimental results, dashed lines to the calculated ones, (the upper pair to v=7 
and the low one to v=8). 
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The angular distributions for the same field characteristics, are gathered in figure ITT1 
They, precisely, correspond to a-dependent ID representations of fixed /c p -cuts of the 3D 
information of the type displayed in figure El This is done for two different values of k p , 
namely k p =5.5 a.u. and k p =6.5 a.u. corresponding to the positions of the maximum am- 
plitudes of the two major peaks of figure HO attributed to v—7 and 8, respectively The 
following observations can be pointed out: 

i) These angular distributions, although labeled as v—7 and 8, actually contain informa- 
tions originating from v =9,10,..., through the red-tail contributions of these levels (as is clear 
from figure El panel c). The excellent theory- versus-experiment agreement that is reached 
has to be judged within this intricate influence of the higher energy part of the spectrum. 
It is also worthwhile noting that due to larger experimental errors affecting v=9,10, angular 
distributions are not studied for higher v 's than 8. 

ii) v =7 is much better aligned than v=8. This is basically due to the bond softening 
mechanism. A high potential barrier at 6 = tt/2 protects v=7 population against pho- 
todissociation. This barrier is lowered at 9=0 or tt where the radiative coupling is at its 
maximum, leading to efficient alignment, that is even better for increasing intensity. In 
other words the wavepacket associated with v=7 has to skirt around a high potential barrier 
at 9 = 7r/2 before dissociating, whereas the one associated with v=8 being closer to the top 
of the barrier can more easily tunnel. The consequence is that dissociation is facilitated for 
8=0, or 7r, when the initial population lies on v=7. 

A better understanding and interpretation of the way following which these complemen- 
tary bound softening and vibrational trapping mechanisms ultimately affect the dissociation 
process, could be gained by a dynamical investigation. Figure ^] illustrates a time-resolved 
decay dynamics of individual vibrational levels. The lower panel gives the temporal shape of 
the laser intensity for the strongest field into consideration (E =0.7 mJ, with the parameters 
of the last row of Table H} . The decay dynamics is given in terms of the decrease of the 
short range population as a function of time, starting from a given v of the parent ion 
(i.e. the time dependence of the norm of the internal region wavepacket ||$ 7 || 2 , as defined 
in subsection 2.1). The behavior of different v's can be classified as follows: 

• Levels affected by the bond softening mechanism; namely v =6,7,8. The decay starts 
only after the maximum of the laser pulse, which is required for the potential barriers to be 
sufficiently lowered. Although the decay mechanism is rather fast (the slope of ||$ / || 2 versus 
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FIG. 12: Time-resolved decay dynamics of individual vibrational levels in correspondance with the temporal 
laser intensity distribution. 
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time is large), the photodissociation starting from v=Q is not complete, due to the fact that 
the potential barrier is closed before total escape towards the asymptotic region. 

• v—9 which lies at the curve crossing region dissociates completely and faster than all 
other levels; 

• Levels affected by the vibrational trapping mechanism; namely v =10,11,12. The popu- 
lations of these levels start to dissociate during the laser rise time, but about the maximum 
intensity their decay rate is lowered (the slope of ||$ 7 || 2 versus time lower than the one of 
levels w=6,7,8). This is basically due to the fact that they are vibrationaly trapped in the 
temporarily closed upper adiabatic potential. It is also interesting to note that the popula- 
tion of v—12, trapped during the radiative interaction, partially returns back to the ground 
state bound potential, in such a way that dissociation starting from v =12 is not complete. 

The dynamical alignment characteristics are illustrated in figure Here again the lower 
panel indicates the temporal shape of the strongest laser. The upper panels display the 
average value of (cos 2 9) on the internal region wavepacket; i.e. ($ 7 | cos 2 0|$ / )/($ / |$ / ), 
indicating the alignment characteristics of the parent ion as a function of time. Three 
initial levels are in consideration, each pertaining to one of the previously selected classes. 
The bond softening mechanism leading to the dissociation of v—6 (panel a) results into 
very efficient alignment during the pulse, which even remains during the fall off regime. A 
thorough interpretation, already given in the literature j^], can be summarized by referring 
to a three dimensional representation of the adiabatic potential energy surfaces (displayed 
in figure EJ) . This is provided in figure El for a single photon dressed ground and excited 
states of including (by diagonalization) the radiative interaction. The photodissociation 
dynamics starting from v—6 is described by a wavepacket that evolves on the lower adiabatic 
potential surface. It has first to skirt around the potential barrier at 9 = n/2 and end up in 
the lower energy valley at 9 = or n. This is why fragments originating from a parent ion 
in an initial state well protected by a hardly penetrable potential barrier (as for v =6,7,8) 
are well aligned through the bond softening mechanism. On the opposite situation are 
parent ions in an initial state basically pertaining to the upper adiabatic potential energy 
surface (as for t>=10,ll,12). This surface presents a minimum around 9 = n/2, where the 
population is temporarily trapped. The dissociation by a single photon absorption further 
proceeds by a non-adiabatic transition to the lower adiabatic surface. Such a transition is 
more efficient for 9 ~ n/2 (induced by a lower radiative coupling). Although the last step 
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FIG. 13: Average values of < cos 2 8 > for the different vibrational levels: (a) v=6; (b) v=9; (c) v=12 ;(d) 
- laser temporal profile. 

which is the evolution on the lower adiabatic surface tends to align the fragments, this effect 
is less efficient as the parent ion is prepared at 9 ~ n/2 on this surface. The result is clearly 
understandable in terms of this vibrational trapping mechanism for v—12 (panel c, figure 
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FIG. 14: 3D adiabatic representation of the potential energy surfaces in the presence of the field with 
peak intensity I = 7.5 x 10 12 W/cm 2 . 

fT3*|) . During the rise time of the laser pulse, the parent ion is well trapped on the upper 
adiabatic surface leading to a misalignment (the bump of the (cos 2 9) curve at about t=550 
fs). There is no noticeable alignment during the whole duration of the pulse. A second 
misalignment, probably due to the nonadiabatic jump, is obtained at t=800 fs. Finally v=9 
which is basically not affected neither by bond softening, nor by vibrational trapping, does 
not show any alignment characteristics as is clear from panel b. 



IV. CONCLUSION 

Once the competition between multiphoton ionization and dissociation processes is dis- 
carded by preparing the parent ion through an electric discharge experiment, a rather 
simple and complete quantum modelisation is provided for a theory-versus-experiment com- 
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parison of angularly resolved kinetic energy spectra of photofragments resulting from intense 
field dissociation. An Abel transformation relates the probability V(k, 6) for a single 
molecule in a given initial ro- vibrational state, to dissociated with kinetic momentum k along 
its polar direction 9 with respect to the laser polarization vector, to P(k p ,a), the proba- 
bility for the photofragment H to be detected on a pixel of the detector plate labeled by 
its polar coordinates (p, a), A quantitative reproduction of experimental data requires some 
statistics over initial ro- vibrational states on one hand and over the spatial distribution of 
laser intensities interacting with molecules positioned at different places in the ionic beam, 
on the other hand. 

An excellent agreement is obtained with experimental spectra and especially for the align- 
ment characteristics of the photofragments. A thorough interpretation can be conducted for 
single vibrational peaks of the spectra in terms of basic mechanisms, such as bond softening 
and vibrational trapping. The most striking observation is the major role the laser volume 
effect is playing, in particular over lower vibrational levels. 

Among our feature prospects, is the elucidation of the role of isotope effects in the pho- 
todissociation of and HD + , which are currently studied in H.Figger's group. 
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Appendix 



This appendix is devoted to some geometrical and and vector relations illustrated in 
figure 121 T*i and r 2 are the vectors pointing H and H + in the laboratory frame, R and Rq 
defined by 

R = r 1 -r 2 (A-l) 



1/ 

Rg = vtuy = -{n +r 2 ) 



(A-2) 



are the relative internuclear separation and the position of the center of mass G, respectively 
(by neglecting the contribution of the electron). The right-angle triangles 00' M and O'M'M 
lead to: 



p 2 + (D- vtf 



and 



|r|| 2 = p 2 + D 2 



whereas from the triangle MGO one gets: 



1 



-R = r-Ra. 



(A-3) 
(A-4) 

(A-5) 



The unitary vector ur along R, can be easily evaluated using Eqs (|A-2IA-3|) and Eq. ()A-5|) : 

R r - (itiy 

UR= mr w+v-vtni* (A " 6) 

and its projection over u y is nothing but: 



U R ■ Uy 



ru„ — vt\\u,. 



D-vt 



(A-7) 



[p 2 + {D - vt) 2 } 1/2 [p 2 + {D - vt) 2 } 1/2 
taking into account: r u y = D. 

The polar angles 9 and a positionning H (and M) in the center of mass and detector 
frames can be related using the right-angle triangle GM'M: 

HH' p cos a 



cos 9 



hR\\ \\^R\ 



or finally taking into account Eq. ljA" 



cos 9 



pcosa 



[p 2 + (D - vt) 2 } 



2ll/2 • 



(A-8) 



(A-9) 
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